NMR relaxation and resistivity from rattling phonons in pyrochlore superconductors 
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We calculate the temperature dependence of NMR relaxation rate and electrical resistivity for 
coupling to a local, strongly anharmonic phonon mode. We argue that the two-phonon Raman 
process is dominating NMR relaxation. Due to the strong anharmonicity of the phonon an unusual 
temperature dependence is found having a low temperature peak and becoming constant towards 
higher temperatures. The electrical resistivity is found to vary like at low temperatures and 
following a \/r behavior at high temperatures. Both results are in qualitative agreement with 
recent observations on /3-pyrochlore oxide superconductors. 
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The recent discovery of superconductivity in the family 
of pyrochlore oxides KOS2O6, RbOs206, and CsOs2 06 
has attracted great interest because of their unusual 
properties. Among these KOS2O6 with the highest Tc 
of 9.6 K appears to be the most unusual. The tempera- 
ture dependence of its electrical resistivity shows a strong 
concave-downward temperature dependence [11,01) in 
contrast to the other two compounds, where a temper- 
ature dependence at low temperatures has been observed 
0) S]- Specific heat measurements have shown a large 
mass enhancement, large specific heat jump at Tc, and 
existence of low frequency Einstein modes [HI, Q . Band- 
structure calculations have indicated that the anomalies, 
in particular in KOS2O6, might be due to a highly anhar- 
monic low frequency rattling motion of the alkali-ion in- 
side an oversized cage formed by the Os and O ions 0, Q ■ 
This is consistent with X-ray observations of anomalously 
large atomic displacements for the K ions P and low 
frequency phonon structures seen in photoemission spec- 
troscopy [io|. 

Recent observations of NMR relaxation rates l/TiT at 
the K site have been demonstrated to be entirely domi- 
nated by the vibrations of the K ion via coupling of the 
electric field gradient to the nuclear quadrupole moment 
[m . Such a domination of phonons for nuclear relaxation 
usually occurs in diamagnetic insulators, but is extremely 
rare in metals and was attributed to the rattling motion 
of the K ions. The temperature dependence of 1/TiT was 
found to be anomalous as well, showing a peak around 
12-14 K and decreasing at higher temperatures. It has 
been argued that such a behavior is inconsistent with 
the two-phonon Raman process, which usually dominates 
quadrupolar relaxation, and it has been interpreted in 
terms of the direct phonon process with a strongly tem- 
perature dependent phonon damping rate [ll[. In the 
superconducting state 1/TiT exhibits a sudden decrease, 
suggesting a strong coupling of the phonon mode to the 
conduction electrons and an associated increase of the 
phonon lifetime in the superconducting state. 



Motivated by these experimental findings, in the 
present work we study the influence of a local, strongly 
anharmonic, and damped phonon mode on the NMR re- 
laxation rate and the electrical resistivity. We find that 
within this model the two-phonon Raman process is ex- 
pected to dominate NMR relaxation over the direct pro- 
cess. Due to the anharmonicity of the phonon mode the 
temperature dependence of 1/TiT for the Raman pro- 
cess is qualitatively different from harmonic phonons. We 
show that the experimental data are qualitatively repro- 
duced, showing a peak at about half the low tempera- 
ture phonon frequency and decreasing towards a constant 
value at higher temperatures. Within the same model 
also the temperature dependence of the resistivity can be 
understood, following a law at low temperatures and 
a -s/T behavior at high temperature. Thus, a qualitative 
understanding of both the anomalous NMR relaxation 
and electrical resistivity can be obtained simultaneously 
by considering the anharmonicity of the phonon involved. 

As a description of the local alkali-ion motion we start 
from an anharmonic Hamiltonian of the form 
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where x, p, and M are the spacial coordinate, momen- 
tum, and mass of the alkali-ion, respectively, and a and 
6 > are constants. We note that according to band- 
structure calculations, in KOS2O6 the quadratic term 
even becomes negative a < 0, resulting in a shallow dou- 
ble well potential, while in the other two compounds a is 
positive [31 ■ We treat this Hamiltonian in a selfconsistent 
quasi-harmonic approximation resulting in an effective 
harmonic Hamiltonian 



H = 



P 

2M 



(2) 



where the effective phonon frequency ujq becomes tem- 
perature dependent and has to be determined selfconsis- 
tently from the equation Mlo^ = a + b{x'^)i^g^T- Here, 
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FIG. 1: Effective phonon frequency uq/ujoo as a function of 
temperature T/uoqo for anharmonicity parameter 13 — 0.1, 1, 
and 10 (double logarithmic scale). 



{x^)u]o,T denotes the thermal average of that depends 
on both Wo and temperature T and is given by 
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Note, that the effective phonon frequency loq is a ther- 
modynamical average frequency and not to be confused 
with the discrete, temperature independent energy levels 
of the Hamiltonian Eq. ([1]). In the zero temperature limit 
we have 
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Note that this equation guarantees ujo{T = 0) > for 
b > 0, even if a < 0. In the high temperature limit 
T ^ Wo we find 
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We can eliminate the parameter a in favor of loqq — 
ujo{T — 0) using Eq. (|4]). Then we can derive the fol- 
lowing nonlinear equation: 
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where /3 = > is a dimensionless parameter char- 



acterizing the amount of anharmonicity. In Fig. [1] is 
shown as a function of T/ujqo for /3 = 0.1, 1, and 10. 
monotonously increasing function of temperature. 
Defining the retarded phonon propagator as 



IS a 



D{lu) = 



2woM 



{[x{t),xm) (7) 




(a) 



(b) 



FIG. 2: Two lowest order processes leading to a finite phonon 
lifetime: (a) coupling of the phonon to the conduction elec- 
trons and (b) fourth order phonon process due to the anhar- 
monicity. 



the interaction of the phonon with the conduction elec- 
trons can be described by the phonon self-energy 11 (w) 
via Dyson's equation D~^{uj) = Dq^{uj) — n(w). Here, 
the noninteracting phonon propagator due to the effec- 
tive Hamiltonian Eq. ^ is given by Dq{ui) = ^f^ji 
Thus, the phonon propagator can be written in terms of 
the self-energy as 



D{^) = 
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The real part of the self-energy leads to a renormaliza- 



tion of the phonon frequency w^ 



2woRe H(w) 



while the imaginary part determines the damping width 
r(w) = — Im H(a'). Due to symmetry the damping width 
has to be an odd function r(w) = — r(— w), while the 
real part has to be even Re H(w) = Re H(— w). In Fig. [2] 
we show the two lowest order processes contributing to a 
finite phonon damping. Process (b) is due to the fourth 
order process neglected in Eq. ([2]). This process does not 
lead to a change of the phonon damping, when supercon- 
ductivity sets in. As mentioned above, the NMR results 
suggest that there is a strong coupling of the phonon to 
the conduction electrons and therefore we expect process 
(a) to play the leading role. Since the Fermi energy of 
the conduction electrons is a much larger energy scale 
than the phonon energies, process (a) will be featureless 
in the low energy range of interest in the normal state. 
The damping width is a linear function at low energy 
and we can write T{oj) — auj. The real part of the self- 
energy can be assumed constant and must be negative 
(i.e. the electron-phonon interaction leads to a softening 
of the phonon). The dimensionless parameter a deter- 
mines the phonon damping and is given by a = g'^N^{0), 
where N{0) is the density of states at the Fermi level 
and g the electron-phonon coupling constant. It is im- 



portant to keep in mind that g 
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because here 



Wo is temperature dependent. Therefore, in the follow- 
ing we write a = Tq/ujq, where Fq is a (temperature 
independent) phonon damping rate and will be treated 
as a phenomenological parameter. The phonon spectral 
function is given by 
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This is equivalent to the famihar damped harmonic oscil- 
lator formula apart from the fact that tOr becomes tem- 
perature dependent here due to the anharmonicity. 

The quadrupolar NMR relaxation rate is usually dom- 
inated by the two phonon Raman process [l2i] . which is 
given by 
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where V2 is proportional to the second spacial derivative 
of the electric field gradient and col is the nuclear Lamor 
frequency, which is usually very small. In Ref. [Til it 
has been suggested, however, that in the case of KOS2O6 
the direct (one phonon) process might be dominating in- 
stead. Here, we are comparing these two processes for 
our model. For the direct process we have 
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where Vi is proportional to the first spacial deriva- 
tive of the electric field gradient. Using the 
fluctuation-dissipation theorem we can directly re- 
late Tjvrj to the phonon spectral function: tttu = 

'^^Vi 2J"oM i^e-'-^^L/'^BT M^L)- Assuming hujL « keT 
we find 
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As has been pointed out in Ref. Illl this expression is 
monotonously decreasing as a function of temperature 
and thus cannot reproduce the observed peak in jr^- 
Due to the temperature dependence of ujq at high tem- 
peratures we have ^tdy ^ '-^q'^ ^ T^^. 

Following the discussion of the Raman process in 
Ref. JJ, we can also express Eq. (flUl) in terms of the 
phonon spectral function: 
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(13) 

At sufficiently high temperatures T ^ LOr, in the normal 
state we can approximate the Rose function n(a;) ~ 
and Eq. can be evaluated analytically. After some 
algebra we arrive at 



1 _ (kBTV2V 2tjg4r2 



\ 2uJoM J To uj$ 

Due to the temperature dependence of u>o in Eq. 
find the high temperature behavior 
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Thus, TfTTiTp approaches a constant high temperature value 
for anharmonic phonons. This is in strong contrast to 
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FIG. 3: Temperature dependence of 1/T\T for the Raman 
process and anharmonicity parameter {3=2 (dotted), 5 
(dashed), and 20 (solid). The dashed-dotted line shows the 
direct process for comparison. 



harmonic phonons, where =^ is increasing linearly with 
T. To estimate the relative importance of the two pro- 
cesses we can use Eqs. (fT4|) and (fT2|) to calculate the ratio 
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From this expression it becomes clear that the Raman 
process is certainly dominating at high temperature and 
also, when the phonon damping Fq is small. We note that 
the ratio ^ has dimension of an inverse length, which 
should be of the order of magnitude of the size of the 
potential, in which the K ion moves. Therefore, as an 

order of magnitude estimate we assume t-%- — ^ 1 

In the normal state at low enough temperatures 
Eq. nil) leads to 
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because of the linear oj dependence of Aiw). As an il- 
lustration Fig. [3] shows a numerical calculation of ^ 

-'i ^ 

as a function of T/luoq based on Eqs. and ^ for 
Tq/ujoo = 0.1, Re n = 0, and three different values for the 
anharmonicity parameter /9 = 2, 5, and 20. Clearly, with 
increasing anharmonicity a low temperature peak devel- 
ops near T w luqq/2 and becomes more pronounced. This 
can be understood from the uJq^ behavior in Eq. (llSp : at 
low temperatures luq becomes constant and all curves fall 
on top of each other. However, at higher temperatures 
ujo increases the stronger the larger (3 is, resulting in a 
quick decrease for larger f3. For comparison, the dashed- 
dotted line shows the direct process for /3 — 20, assuming 
1. Clearly, the direct process can be ne- 
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FIG. 4: Temperature dependence of the resistivity for the 
same parameters as in Fig. [31 Inset: Temperature dependence 
of 1/T\T fitted to the experimental data of Ref. [ll] (sample 
A) on a double logarithmic scale. 



that the direct process is exponentially suppressed in the 
superconducting state. 

In the inset of Fig. 2] we show a fitting of Eq. to 
the experimental data of Ref. [ll| (sample A) . From this fit 
the following parameters are found: ujr{T = 0) = 27.4 K, 
/3 = 6.27, Fo = 4.0 K, and u^Re H =-(21 Kf . The low 
temperature value of the renormalized phonon frequency 
bJr{T = 0) = 27.4 K is in reasonable agreement with a 
low frequency Einstein frequency of 22 K found in fits to 
specific heat data [6] and 24.4 K found from a phonon 
structure in photoemission 10]. 

In the following we calculate the resistivity due to cou- 
pling to the local anharmonic phonon Eq. We will 
follow the work by Mahan and Sofo [l3| . The imaginary 
part of the retarded electronic self energy is given by 
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is the Fermi function. The re- 
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where 



where f{fl) 

sistivity can be calculated from p{T) 
Ljp is the Plasma frequency and the electronic lifetime is 
given by t{T) = - J^^dE t{E)-^. We note that we 
can neglect vertex corrections here, because the phonon 
is local, i.e. momentum independent. Therefore, the 
transport lifetime is the same as the electronic lifetime. 

Because A{il.) is linear in at low frequencies, we find 
from Eq. (fTS)) that at temperatures T uir the resistivity 
varies like p T"^. However, in the high temperature 
limit we find from Eq. (fH 
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Therefore, at high temperatures the resistivity varies as 
p ~ Vt, a concave-downward temperature dependence. 
In Fig. [4] numerical results for the resistivity p are shown 
for the same parameters as in Fig. [3l With increasing /3 
the concave-downward temperature dependence becomes 
more pronounced. Considering woo ^ 25 K this result 
is in qualitative agreement with the resistivity data on 
KOS2O6 [1, 0]. A quantitative fitting would require tak- 
ing into account the other, harmonic phonon modes as 
well. Since at present not much is known about these a 
more quantitative comparison has to await more detailed 
experimental information. 

In conclusion, we have calculated the temperature de- 
pendence of NMR relaxation rates and resistivity due to 
coupling to a strongly anharmonic phonon mode. Our 
results are in agreement with recent experiments on /3- 
pyrochlore superconductors, where such an anharmonic 
rattling motion of the alkali-ions is believed to cause 
various anomalies. Our theory may also be applica- 
ble to other rattling systems like the filled skutterudite 
LaOs4Sbi2, as a recent NMR study suggests [14] . 

Thanks are due to M. Takigawa and Z. Hiroi for valu- 
able discussions and for providing their data. 
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